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1.  Justification  of  the  basis  for  weighted  rankings  analysis 


The  method  of  weighted  rankings  analysis  was  introduced  in  Quade  (1972, 

1979).  A  brief  description  is  as  follows:  Let  be  the  yield  of  the  j-th  treatment  in 

the  i— th  block,  for  i  =  l,...,n  and  j  =  l,...,m  and  let  Rjj  be  the  within— block  rank  of 

Xjj.  Let  Dj  be  a  measure  of  apparent  variability  for  the  i-th  block,  Q.  be  the 

corresponding  rank.  Let  t^,...,tm  be  any  constants,  and  let  Sp...,sn  be  constants  such 

that  0  <  st  <  ...  <  s  .  Then  a  statistic  of  the  form 
-  1  -  -  n 


c„  = 


m  f  n  w 

S  S  S0  *R 
j=l  4=1  Qi  RijJ 
n  0  m  0 
53  S.  53  t. 
i=l  1  i=l  J 


may  be  used  for  testing  the  hypothesis  of  no  treatment  effects. 

The  method  of  weighted  rankings  is  based  on  the  idea  that  if  some  blocks 
appear  more  variable  than  others,  then  they  are  perhaps  better  referred  to  as  more 
discriminating.  Hence  it  seems  intuitively  reasonable  that  these  blocks  receive 
greater  weight  in  the  analysis. 

In  this  research  we  examine  this  credability  hypothesis  and  provide  objective 
evidence  showing  how  the  more  variable  blocks  do  indeed  better  reflect  any  true 
treatment  effects.  The  n— block  two  treatment  exponential  case  was  considered,  and 
the  foundation  is  established  to  extend  the  result  to  the  more  general  cases  (see 
Appendix  I). 

2.  A  nonparametric  multivariate  test  for  homogeneity  based  on  all  nearest  neighbors 

Let  {X.j:  i  =  1,2;  j  =  l,...,n.}  be  two  independent  random  samples  of 
observations  in  the  Euclidean  space  R^,  where  X-j  has  distribution  function  Fj.  The 
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problem  under  consideration  is  to  test  the  hypothesis  Hq  :  =  F2,  against  the 

completely  general  alternative  Ha  :  Fj  f  Fj. 

Let  ROJji'J')  be  the  rank  of  observation  X. ,  j,  with  respect  to  nearness  to 
X..;  we  assume  that  there  will  be  no  ties.  Then  we  define  X.,j,  as  "the"  k— th  nearest 
neighbor  of  Xy  if  R(i,j;i/ ,j' )  =  k,  and  as  "a"  k-nearest  neighbor  if  R(i,j;i'  ,j' )  <  k. 
Interest  in  statistical  procedures  based  on  such  nearest  neighbors  has  grown  as 
high-speed  computers  have  made  the  application  of  these  techniques  practicable, 
since  the  idea  of  making  inferences  about  an  object  based  on  nearby  objects  appears 
to  be  a  fundamental  mechanism  of  human  perception.  A  review  of  early  work  using 
nearest— neighbor  approaches  to  our  problem  may  be  found  in  Schilling  (1986). 

Schilling’s  own  approach  is  as  follows.  Let 


Lj(r)  =  I{r-th  nearest  neighbor  of  X.j  is  in  sample  i} 
for  i  =  1,  2;  j  =  1 . n^ 


r  =  1,...,N— 1  (N  =  n^+ng), 

where  I{E}  is  the  indicator  function  of  the  event  E.  Count  the  number  of  k— nearest 
neighbors  to  X-j  which  are  in  the  i-th  (same)  sample,  viz. 


c*  “  l 


Summing  these  counts  over  all  observations  yields  what  may  be  called  the  "Schilling 
total",  of  order  k: 


His  test  statistic  for  our  problem  is  then 


Sk  =  Tk/Nk, 

which  is  "the  proportion  of  all  k-nearest  neighbor  comparisons  in  which  a  point  and 
its  neighbor  are  members  of  the  same  sample"  (1986,  p.  800),  or  "Schilling 
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proportion"  of  order  k.  One  would  expect  to  have  a  larger  value  under  H&  than 
under  Hq  because  of  a  lack  of  complete  mixing  of  the  two  samples  when  the  parent 
distributions  are  not  identical;  hence  large  values  of  may  be  considered  significant. 
Schilling  shows  that  the  asymptotic  distribution  of  under  Hq  is  normal,  for  every 
positive  integer  k,  and  that  the  test  which  rejects  for  large  values  of  is  consistent 
against  the  general  alternative  H&. 

Schilling’s  work  suggests  that  the  choice  of  order  is  not  of  great  importance; 
nevertheless,  it  is  arbitrary,  and  he  gives  no  guidance  for  choosing  it.  To  resolve  this 
issue  we  propose  to  take  as  test  statistic  the  sum  of  the  Schilling  totals,  which  is 
equivalent  to  a  certain  weighted  average  of  the  Schilling  proportions: 

W  =  ETk  =  NIkSk. 

This  statistic  may  have  intuitive  appeal  in  that  it  is  equivalent  to  a  sum  of  N 
Wilcoxon  rank  sums,  and  is  also  a  linear  combination  of  two  U-statistics,  as  we  have 
shown.  We  study  some  of  the  exact  properties  of  W.  For  the  asymptotic  properties, 
we  note  that  W  is  equivalent  to  a  sum  of  Schilling  proportions,  each  of  which  is 
asymptotically  normal;  but  asymptotic  normality  of  their  joint  distribution  has  not 
been  shown.  Similarly,  we  have  shown  that  W  is  a  sum  of  Wilcoxon  statistics,  each 
of  which  is  asymptotically  normal;  but  again  asymptotic  normality  of  their  joint 
distribution  has  not  been  shown.  We  might  also  attempt  to  show  this  by  means  of 
U-statistics  theory,  but  have  not  yet  been  able  to  work  out  the  necessary  conditions 
on  the  variance.  Nevertheless,  it  seems  reasonable  to  conjecture  that  W  itself  is 
asymptotically  normal.  Computer  simulation  was  used  to  estimate  the  95th 
percentile  of  W,  then  this  value  will  be  used  as  the  critical  value  for  a  test.  That  is, 
the  null  hypothesis  will  be  rejected  if  W  is  greater  than  this  critical  value.  The 
empirical  power  of  the  test  based  on  W  was  also  calculated  based  on  this  critical 
value. 
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3.  Exact  tables  for  Spearman’s  footrule 

Spearman  (1904)  based  the  measure  of  rank  correlation  as  his  "footrule"  on 

D(<7,*)=  E  |o(i)-*(i)|, 
i=l 

where  a  and  r  are  any  two  permutations  of  SQ,  the  set  of  all  permutations  of  the  first 
n  integers.  Ury  and  Kleinecke  (1979)  provided  the  exact  commulative  distribution  of 
D  for  n  =  2(1)10.  Franklin  (1988)  has  extended  the  tabulation  to  n  =  11(1)18.  The 
difficulty  in  extending  these  tables  is  due  to  the  exponential  growth  rate  of  time 
needed  to  generate  all  permutations  of  Sn-  By  relating  D  to  a  Markov  chain  on  Sn, 
we  extend  the  exact  tabulation  to  n  =  19(1)40.  We  also  investigate  the  adequacy  of 
approximation  to  the  normal  distribution. 

4.  The  asymptotic  permutational  normality  of  certain  weighted  measures  of 

correlation 

Let  be  given  observation  (X.,Yj),  i  =  l,...,n,  on  the  continuous  bivariate 
random  variable  (X,Y).  Without  loss  of  generality  we  may  relabel  the  X’s  to  have 
ranks  l,...,n:  then  let  the  corresponding  ranks  of  the  Y’s  be  Rp"-iRn-  Define 

t=  =  i  i{Ri  <  j} 

J  i=i 

where  I{E)  is  the  indicator  function  of  the  event  E,  and  let 

n 

T  =  E  w.T. 
j=i  1  1 

where  w^,...,wn  is  a  sequence  of  weights.  Then  T  may  be  regarded  as  a  weighted 
measure  of  rank  correlation  between  the  X’s  and  Y’s.  This  notion  was  introduced  in 
Salama  and  Quade  (1982),  although  they  discussed  only  the  special  case  where  Wj  = 
1/1 
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In  this  research,  a  martingale  approach  is  used  to  establish  the  asymptotic 
normality  of  a  class  of  weighted  correlation  statistics. 

5.  Spherical  uniformity  and  the  Cauchy  distribution 

To  motivate  the  characterization  problems  (to  be  posed),  we  consider  the 
following.  Let  Y  be  a  stochastic  (m-)  vector  for  which  EY  =  X/?,  where  X  is  a  known 
m*n  (design)  matrix  and  /?  is  an  n— vector  of  unknown  regression  parameters.  Then 

A  A 

P ,  the  least  squares  estimator  (LSE)  of  P ,  satisfies  the  equation:  (X'X)/?,  =  X' Y. 
When  X'X  is  not  of  full  rank,  additional  constraints  are  imposed  on  P  in  order  that 
the  LSE  may  be  defined  uniquely.  In  the  absence  of  a  physically  natural  form  of  such 
a  (linear)  constraint,  often  this  choice  is  made  rather  arbitrarily.  In  some  random 
effects  or  mixed  models,  when  Rank(X'X)  =  n— 1,  an  additional  random  constraint  is 
taken  as  c '/?  =  0,  where  ||c|j  =  1.  In  robust  regression  (cf.  Huber  (1981,  p.  170)), 
when  n  is  large,  often  c  is  chosen  at  random  with  respect  to  the  invariant  measure  on 
the  unit  sphere  ||c||  =  1.  In  the  context  of  asymptotic  normality  of  estimates,  Huber 
(1981)  assumed  that  Rank(X'X)  =  n,  and  we  shall  see  that  a  different  picture 
emerges  when  Rank(X'X)  is  less  than  n.  Suppose  that  Rank(X'X)  =  n-1  and  let 
c'/?  =  0  be  an  additional  constraint  such  that  ||c[|  =  1.  Let  /?  be  the  solution 

V 

A 

satisfying  (X'X)/?  =  X' Y  and  c'/?  =  0.  Let  v  be  the  eigenvector  corresponding  to 

A 

the  null  eigenvalue  of  X'X.  Then  /?  may  be  wriaen  as 

V 

\  =  \  +  ((c/^oMc/v))v’ 

A 

where  /?q  is  any  particular  solution  of  (X'X)/?  =  X' Y  (and  is  therefore  a  random 

A 

variable).  Thus,  given  /?q,  we  may  write 

(0c-0q)  =  vl(c'0o)/(c'v)],  (2) 

A 

and  its  conditional  distribution,  given  /?q,  is  generated  exclusively  by  the  uniform 
distribution  of  c  on  the  spherical  surface.  Note  that  v  is  nonstochastic  (as  X'X  is),  so 
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that  the  key  factor  on  the  right-hand  side  of  (2)  is  the  ratio  (c'/?q)/(c'v),  where  0q  is 
held  fixed.  This  leads  us  to  the  following  problem. 

For  two  arbitrary  points  b  and  v  in  Rn,  let  L(t)  =  b  +  tv,  t  €  R  be  a  given 
line.  Let  Sn_1  =  {c:c'c  =  1}  be  the  unit  sphere,  and  for  every  c  6  Sn_1,  let  P(c)  be 
the  hyperplane  defined  by  P  c  =  0,  P  6  Rn.  Let  L(t  )  =  b  +  (t  )v  be  the  point  at 

V  V 

which  L  and  P(c)  intersect,  for  c  6  Sn— .  Define  then 

=  x„(c) =  sten(tc)llL(tc) —  Mo)ll  —  8>6n(tc) I *c I IMI.  ces”-1.  (3) 

Assuming  that  c  has  a  uniform  distribution  on  the  sphere  Sn— \  we  show  that  for 
every  n  >  2,  Xn  has  the  same  distribution  as  aX  +  d,  where  a  and  d  are  real  numbers 
and  X  has  the  standard  Cauchy  distribution.  Some  other  related  characterization 
results  are  also  considered  in  the  same  view. 

6.  Topological  entropy  of  countable  Markov  chains 

We  consider  a  symbolic  dynamical  system  (X,a)  on  a  countable  state  space. 
We  introduce  a  kind  of  topological  entropy  for  such  systems,  denoted  h*,  which 
coincides  with  usual  topological  entropy  when  X  is  compact.  We  use  a  pictorial 
approach,  to  classify  a  graph  T  (or  a  chain)  as  transient,  null  recurrent,  or  positive 
recurrent.  We  show  that  given  0  <  a  <  0  <  a>,  ihere  is  a  chain  whose  h*  entropy  is  0 
and  where  Gurevic  entropy  is  a.  We  compute  the  topological  entropies  of  some 
classes  of  chains,  including  larger  chains  built  up  from  smaller  ones  by  a  new 
operation  which  we  call  the  Cartesian  sum. 
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APPENDIX  I 


REMARKS  ON  JUSTIFYING  THE  INTUITIVE  BASIS 
FOR  THE  METHOD  OF  WEIGHTED  RANKINGS 

Ibrahim  A.  Salama,  School  of  Business,  North  Carolina  Central  University,  Durham 
Dana  Quade,  Department  of  Biostatistics,  University  of  North  Carolina  at  Chapel  Hill 


1.  INTRODUCTION 

Let  be  the  observation  of  the  j-th  of  m 
treatments  in  the  i-th  of  n  complete  blocks,  and 
consider  the  hypothesis  of  no  treatment  effects, 
specifically 

H0:  X  n,  . . .,  X jm  are  interchangeable  for  each  i. 

(By  definition  random  variables  are 
interchangeable  if  their  joint  distribution  function 
is  invariant  under  permutation;  this  implies  that 
they  have  identical  marginal  distributions  and 
equal  —  but  not  necessarily  zero  —  correlations.) 
We  assume  throughout: 

(I)  Independent  Blocks 

For  i=l,  n,  the  random  vectors  X,=(X,,, 

. . .,  X,m/  (the  blocks),  are  mutually 
independent. 

To  simplify  the  exposition  it  is  convenient  also 
to  assume: 

(II)  No  Within-Blocks  Ties 
P{X,V=X  ,}=0forj*j' 

•j 

Thus  with  probability  1  there  will  be  no  ties 
within  blocks. 

The  alternatives  under  consideration  can  be 
fairly  general:  however,  we  have  particularly  in 
mind  that  there  may  be  additive  treatment 
effects,  as  follows: 

UNORDERED  CASE 

H«(u]:  There  exist  quantities  r,,  ...,  rm 
(treatment  effects)  not  all  equal  to  zero, 
such  that  for  i  =  l,  ...,  n,  Xn  —  r,,  ..., 
X/m  —  rm  are  interchangeable. 
ORDERED  CASE 

H«[o]:  The  quantities  r,,  ...,  rm  (as  above) 
satisfy  r,<...<rm,  with  rtjhrm- 

Standard  nonparametric  procedures  for 
attacking  this  problem  are  based  on  within-block 


rankings:  for  example,  the  tests  of  Friedman 
(1937)  and  Brown  and  Mood  (1951)  for  Ha[u]; 
and  Lyerly  (1952),  Page  (1963)  and  Jonckheere 
(1954)  for  lIa[o].  Let  C  ,  be  some  measure  of 
rank  correlation  between" block  i  and  block  i'. 
Then 

e  =  s  ,c.,/(S) 

«/  ••  z 

is  the  average  internal  rank  correlation,  and  we 
reject  H0  in  favor  of  Ha[u]  if  C  is  too  big.  (This 
is  equivalent  to  Friedman's  test  if  Spearman 
correlation  is  used.)  Similarly  let  C{  be  the  rank 
correlation  between  block  i  and  the  ordering 
given  by  the  alternative.  Then 

C  =  ?  CJn 

is  the  average  external  rank  correlation,  and  we 
reject  H0  in  favor  of  Ha[o]  if  C  is  too  big.  (This 
is  equivalent  to  Page’s  test  if  Spearman 
correlation  is  used.) 

The  standard  procedure  is  parametric  two- 
way  analysis  of  variance,  but  this  adds  two 
assumptions: 

(III)  Additive  Block  Effects 

There  exist  quantities  /?,,  ...,  /?„  (block 
effects)  such  that  the  random  vectors 
(Xj,-/?,,  X,m-/3,)'  are  identically 

distributed. 

(IV)  Normality 

The  X^-’s  are  [jointly]  normal. 

By  Assumption  III,  comparisons  of  observations 
are  possible  between  blocks  as  well  as  within,  so 
procedures  which  use  only  within-block 
comparisons  waste  information.  A  method  of 
weighted  within-block  rankings,  which  makes  use 
of  Assumption  (III)  without  requiring 
Assumption  (IV),  has  been  introduced  by  Quade 
(1972,  1979).  The  idea  behind  this  method  is 
that  blocks  in  which  the  observations  are  more 
distinct  are  more  likely  to  reflect  any  underlying 


true  ordering  of  the  treatment  effects.  These 
blocks,  which  may  be  referred  to  as  more 
credible,  should  receive  greater  weight  in  the 
analysis.  (In  practice,  credibility  is  measured  by 
apparent  variability;  but  note  that  by 
Assumption  III  the  true  variability  is  the  same  in 
all  blocks.) 

To  determine  the  weight  for  the  i-th  block, 
use  some  location-free  statistic  D,=D(Xil,  ..., 
Xim)  which  measures  the  credibility  of  the  block 
with  respect  to  treatment  ordering,  and  let  Q,  be 
the  rank  of  D(  among  Dt,  Dn.  Again  for 
simplicity  of  exposition,  make  the  (unessential) 
assumption: 

(V)  No  Between-Block  Ties 

P{D,  =  D  ,}  =  0  for  i  #  i' , 

« 

This  assures  that  there  will  be  no  ties  in  the 
ranking  of  the  blocks.  Let  0<b1<...<b„,  with 
0^bn,  be  a  fixed  set  of  block  scores;  and  weight 
the  i-th  block  proportionally  to  b^  . 

Then  in  testing  against  the  unordered 
alternative  we  use  the  weighted  average  internal 
rank  correlation 

W  =  EE(b0  b0  C  fl/  ((Eb.l’-Eb’) . 

Against  the  ordered  alternative  we  use  the 
weighted  average  external  rank  correlation 

W  =  EbnC,/  Eb,  , 
i=i  Qi  7  <=i  * 

The  purpose  of  this  paper  is  to  examine  the 
notion  on  which  these  weighted  rank  correlation 
coefficients  are  based,  which  we  call  the 
credibility  hypothesis. 

2.  THE  CASE  OF  TWO  TREATMENTS 

With  m=2  treatments,  suppose  the  true 
ranking  is  (1,  2):  i.e.,  r,<r2.  Let  R,^  be  the 
rank  of  X4^  within  the  i-th  block,  and  consider 

P{(R«,  Ria)  =  (1.  2)|Q.=k}  =  (say). 

The  intuitive  notion  is  that 

^  $2  S  •••  ^  * 

where  of  course 


E^j/n  =  <t>  (say) 

is  the  unconditional  probability  P{(Rjl,  Ri2)  = 
(1,  2)).  With  two  treatments  we  may  let  D4  = 
XM-XM.  Then 

(R4j,  R,j)  =  (1,  2)  Dj  >  0 

and 

Qi  =  k  «  |Df |  has  rank  k  among  |D,|,  ...,  |Dn|. 

In  the  very  special  case  of  2  treatments  and  2 
blocks,  let  us  define 

*  =  P{D,  >  0} 

and 

9  =  P{D,+D2  >  0). 

Then 

d>t=  P{D1>0||D1|<|D2|}  =4-(8-4)=24-9 

<h  =  P{Dj>0||D1|>|D2|}  =4>+(0-4>)  =  6 

and  the  credibility  hypothesis  holds  if  0><f>. 
Suppose  (X,j  — r,,  Xi2  — r2)  are  IID  normal  with 
variance  <r3]  then  it  is  easily  seen  that 

where  4  is  the  standard  normal  distribution 
function,  and  0>4>  if  r j<r2.  On  the  other 
hand,  suppose  (X41  —  tx,  Xj2  — r2)  are  IID 
Cauchy;  then  8  =  d>,  and  the  credibility 
hypothesis  docs  not  hold.  Note,  by  the  way, 
that  with  only  two  treatments  the  weighted 
rankings  procedures  are  equivalent  to  signed- 
rank  procedures,  which  are  thoroughly  discussed 
in  Chapter  3  of  Pratt  and  Gibbons  (1981). 

If  we  do  not  limit  ourselves  to  additive 
treatment  effects,  we  may  consider  the 
interesting  special  case  where  X  and  Y  are 
exponentially  distributed,  with  parameters  (say) 
A  and  p,  respectively.  Then  it  is  easily  shown 
that  ^=A/(A  +  /<),  0  =  <j>7(Z  —  24>),  and  0><k  if 
and  only  if  A<p. 

3.  THE  GENERAL  CASE 

We  now  turn  our  attention  to  the  general  case 
where  we  have  m  treatments  and  n  blocks. 


Write 


(R4i» . ..,  R,m)  —  R,  i 

and  consider 

P{R4=£lQ4=k} 

=  nP{R<=|,  Q(=k}  since  P{Qi=k}  =  l/n 

=  W<  D-  <  D'.k} 

=  nPtR^r,  <  D^tj} 

where  D.x  is  the  j-th  order  statistic  from  a 
sample  of  (n— 1)  values  of  D  [that  is,  all  values 
except  D,]. 

Theorem: 

Let  gk_lk  be  the  joint  density  function  of 
and  D^j.  Then 

P{R4=£lQ4=k) 

oo» 

=  n /  /P{R,=r,  a<D,<b}gt_1  t(a,  b)dadb. 
o  o  ' 

Proof: 

P{R.=f|Q4=k} 

=  nP{R,=r,  Q,=k} 

=  nP|R,=r, 

D,j<  •••  <  D,t  i<  D4  <  Dit<  ...  <  Din  l} 

where  (ilt  in_j)  is  a  permutation  of  (1, 
i— 1,  i+1,  n).  This  in  turn  is  equal  to 

nP{R<=£»  ^(*-1)  <  D;  <  d^tj} 

where  d, ..  is  the  j-th  order  statistic  of  a  sample 
(n— 1)  oDscrvations  on  D,  and  thence  equal  to 
the  integral  of  the  Theorem.  QED 

This  form  allows  the  possibility  of  actual 
computation.  The  density  g  is  well-known  given 
the  distribution  of  D.  And  the  probability 
expression  in  the  integrand  depends  only  on  the 
m-dimensional  joint  distribution  of  the 
observations  within  a  single  block.  However,  it 
does  involve  working  with  distributions  of  order 
statistics  from  heterogeneous  distributions. 

In  the  special  cases  of  k=l  and  k=n  we  have 


=  n?  P{R4  =  r,  0  <  Dj  <  t}f(1)(t)dt 

and 

P{R*  =  £|Q.  =  n} 

=  nJPtRj  =  r,  t  <  D<  <  oo)  f(B_1}(t)dt , 

where  f^  is  the  density  function  of  D^. 

4.  AN  EXPONENTIAL  SPECIAL  CASE 

In  this  section  we  consider  an  application  of 
the  preceding  theorem  to  a  simple  and 
analytically  tractable  case.  We  consider  two 
treatments  realized  by  the  random  variables  X 
and  Y,  which  we  assume  independent  with 
probability,  density  functions  f*(x)=e‘*, 
f,(y)  =  Ae  *,  and  0  <  X,  Y  <  oo.  For  the  n 
blocks  we  have  independent  observations 
?4=(Xit  Y,),  i=l,  ...,  n.  Let  D,  =  |X,-Y,.|, 
i=l,  ...,  n.  Assuming  that  A>1,  the  “correct" 
ordering  of  (X,  Y)  is  given  by  the  permutation 
<r=(2,  1).  [The  “correct"  ordering  means 
P(X>Y)  >  P(Y>X)].  (In  this  notation  the 
range  of  the  observations  in  the  i-th  block  is  less 
than  that  of  the  (i+1)*'  block.) 


Since  Corr(R(,),  <r)  =  | 

1  ,f  R(0,i 

=  2 

-1  if  R(i)ii  = 

=  1 

it  follows  that 

E(Corr(R(j),  <r))=  P{R(|))1  = 

2>  P{R(<),i 

=1} 

=  p{R(4),.=2}- 

[i-p{R(.),i= 

=2}] 

=  2P{R(0,,=2}-1  . 

Thus,  we  compute 

Pk  =  P(R(X<)=2|R(D,)=k),  k  =  l,  ...,  n  +  1 

and  show  that  (Pt)  is  strictly  monotone 
increasing  in  k.  The  following  lemma  will  be 
used  in  showing  the  result. 

Lemma: 

Let  X(,j,  ...,  X(nj  be  the  order  statistics  for  a 
sample  of  size  n  from  11  (x),  0  <  x  <1.  Let 
dk=E(X(t+,s— E(X(tj).  If  h(x)  is  monotone 
increasing  (decreasing),  then  (d4)  is  monotone 
decreasing  (increasing)  in  k. 


P{R<)  =  £|Q4  =  1} 


Eiggfr 


similar. 


E(X(t)) 

-H(x)]n‘*dx 

-TJ*(?)H<(x)Il~H(x),n’<dx 

=  1  -  }.E  (j)H<(x)[l-H(x)]n-dx, 

So 


tr,H*-1(x)h(x)[l 

w^xr-xis 


Theorem: 

Let  (x,,  yj),  (xn,  yn)  be  the  observations 
corresponding  to  a  design  with  two  treatments 
and  n  blocks.  Assume  that  X  and  Y  are 
independent,  with  density  functions  f*(x)  =  e'r 
and  f*(y)  =  Ae"  ,  0  <  X,Y  <oo.  Let  Pt  = 
PfRfXJ  =  2|R(D.)  =  k},  where  D,  =  |X,-Y,|. 
If  A>1,  then  {Pt}  is  monotone  in  k;  that  is, 
P,<Pj<...<P«. 

Proof: 

Let  D(X,  Y)  =  |X-Y|,  then 


d*  =  E(X(Jk+1)— E(X(t)) 

=  }(£)H‘(x)[l-H(x)r‘dx. 

Let  y=II(x),  so  x=H*1(y),  and  set  p(y)=dx/dy. 
Then 

d*  =  ^(y)(E)y‘(1-y)B'‘dy 

=/p(y)Lfc(y)dy  (say). 

0 


G(l)=P(D<l)=l-f^,"-I^«-1', 

g(t)=rfA(e',+c'A<)’  0<t<o°- 
We  also  have 

P{R(X)=2,  0<D<t}=r^(l-e-»), 
P{R(X)=2,  t1<D<t2)=I^(e“1-e“2), 

and 

P{R(X)=2,  t<D<oo}=r^e-‘. 


Clearly,  L4(0)=Li(1)=0,  Lt  is  unimodal,  and 
there  exists  a  y*  such  that  0<y*<l  with 
bt(y*)=L  t+1(y*).  Now  if  h(x)  if  monotone 
increasing,  then  /i(y)  is  monotone  decreasing  and 

d*+x=//(y)[Lt+i(y)-L*(y)]dy 

-7  p(y)(Lt(y)-i-*+i(y)]dy 

0 

<  p(y*)/{L*+i(y)-Lt(y)dy 

-  p(y*)7  (L*(y)  -  Li+i(y)]dy 

0 

=  p(y*)/[Lt+1(y)-L*(y)]dy 
0 

=  o. 

(Note  /Lt+l(y)dy=/Lt(y)dy  =  -U  .] 

0  0  *■  i  * 


By  the  preceding  theorem,  for  k=2,  ...,  n  we 
have 

Pt  =  P{R(X)=2|R(D)=k) 

=n°J  / P{R(X)=2,  t!<D<t2} 
o  o 

8(t-l,t)(t2*tl)dtld*’ 

=  n  /  /(f^j)(e  1-e  tl)dtIdt2 

-  "  °?  'ru*'  r'*2!  (n'1)! 

~i+5i  S'*  *  ;(k-2)!(n-l-k)! 

[G(T1»*_2g(t1)g(t  2)[1 — G(t2)]n_1_*dtldtJ 


Hence  {dk}  in  monotone  increasing  in  k.  The  /  «g(*j(t)dt= 
case  where  h(x)  is  monotone  decreasing  is  0 


op-H  (n-1)!  r  A  e-«  A 
'£  (k-l)!(n-l-k)!L  1+5*  1+5*  J 

<i+i>',+e‘ a1iV+i+j'a']  ' 

(°-l)t  n  A  V  1  Xxfl 

i  A(k-l)!(n-l-k)!|_1  1+5  A  1+5  A  J 

*(l+5)(X+xA)[iisX+I+5xA]  (x)1*’ 

(via  the  transformation  X=e'<) 

•/X  Pi  -5  x 

'£  (k-l)!(n-l-k)!|_  I+5A  I+5A  J 

(r&)(l+xAltmx+i+xxA}h‘ 


where  H(x)  is  the  distribution  function  given  by 

h(x)=I+5x+I+5xA’  °^x^1- 

Thus  we  have 

/  «  <8(t)(t)dt=E(X(n_t)), 

where  X ^  is  the  k-th  order  statistic  of  a  sample 
of  size  n-1  from  H(x).  Now 

p‘=(-TTr[E<x(«-*.n-E(x, 

and  for  A>1,  h(x)  is  strictly  monotone 
increasing,  hence  P«>Pn.i  ...  >  Pt. 
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